Calibration of offset, gain and phase errors in m-channel time-interleaved analog-to-digital converters

ABSTRACT

Techniques for correcting component mismatches in an M-channel time-interleaved Analog to Digital Converter (ADC). In order to obtain an error measure for offset, gain or phase, errors, outputs from each ADC are either summed or averaged over N o  samples. Calling each of the sums or averages as X k  where k=1, 2, . . . , M, there are M such values as a result. A single value representing the mean of these M values, X mean , is chosen as a reference value. The offset, gain and phase errors for the M different ADCs are then obtained from X k −X mean . The sign of each offset error, i.e., sign(X k −X mean ), is then used to drive an adaptive algorithm whose output represents an offset correction value for the corresponding ADC. The offset, gain, and phase correction outputs from the adaptive algorithm is fed to an array of Digital-to-Analog converters (DACs) whose outputs are voltages or currents that directly or indirectly controls the offset, gain or phase setting of each individual ADC. Thus, there are M different offset, gain and phase error signals and M different adaptive algorithms operating in conjunction with M different DACs providing offset control signals to M different ADCs. In certain embodiments, spur frequencies can be reduced with the use of notch filters.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 61/233,571, filed on Aug. 13, 2009, and is a continuation of U.S. patent application Ser. No. 12/419,599, filed Apr. 7, 2009, the entire contents of each are also incorporated by reference herein in their entirety.

BACKGROUND OF THE INVENTION

An efficient way of providing very high sample rates, rates that cannot be provided by a single Analog-to-Digital Converter (ADC), is to use a parallel connection of slower ADCs operating in a time-interleaved fashion. An M-channel time-interleaved ADC (MCTIADC) comprises of M ADCs, each operating at a sample rate that is 1/M of the overall system sample rate. In the absence of any impairments or mismatch errors between the ADCs, i.e., assuming all the ADCs are either ideal or have exactly the same characteristics, the output samples appear at equally spaced intervals in a manner that creates a seamless image of a single ADC operating at the system sample frequency.

In practice, however, there are component mismatches between the different ADCs that severely degrade the performance of the MCTIADC system. The commonly offsets and gains of all the ADCs are not the same and the ADCs do not sample at uniform sample instants of the system sample frequency. These mismatches give rise to unnecessary frequency tones or spurs in the spectrum of the signal that significantly reduce the performance of the MCTIADC system. A typical variation of Signal-to-Noise ratio (SNR) is shown in FIG. 1 wherein a tone is swept from a low frequency to almost half the sample rate of the MCTIADC system for various mismatch errors. As can be seen from the figure, the performance of the four-channel ADC is severely hampered due to these errors. Hence, it becomes imperative to estimate and correct these errors to improve the performance of the MCTIADC system.

Herein are shown techniques to minimize the effects of offset, gain and sample-time mismatches by appropriately estimating and correcting these errors in an adaptive manner. In addition, also shown is that the adaptive method can be used in a blind mode wherein the use of any particular calibration signal is circumvented. In other words, the input signal itself serves as the calibrating signal to estimate and correct the mismatch errors.

SUMMARY OF THE INVENTION

This invention generally deals with the estimation and correction of offset, gain and timing errors in an M-channel time-interleaved Analog-to-Digital Converter (MCTIADC). The offset errors which manifest due the difference in errors between the individual ADCs produce spurious frequency content in the spectrum, are called the offset spurs. Assuming that F₅ is the sampling frequency of the MCTIADC system, each ADC samples at the rate of F₅/M and the offset spurs are produced at kF₅/M frequencies irrespective of the frequency or amplitude of the input signal, as can be seen in FIG. 3 for a four-channel time-interleaved ADC. Stated in different words, the offset spurs appear at multiples of the sampling frequency of any single ADC. In order to obtain an error measure for the offset errors, output from each ADC is either summed or averaged over N_(o) samples. Call each of the sum or average as X_(k) where k=1, 2, . . . , M. As can be noticed, there are M such values as a result of the summing or averaging operation over the M channels. A single value representing the mean of these M values, say X_(mean), is chosen as a reference offset value. The offset errors for the M different ADCs are obtained as X_(k)−X_(mean). The sign of each offset error, i.e., sign (X_(k)−X_(mean)), is used to drive an adaptive algorithm whose output represents an offset correction value for the corresponding ADC. The output from the adaptive algorithm is fed to a Digital-to-Analog converter (DAC) whose output is a voltage or current that directly or indirectly controls the offset setting of each ADC. Thus, there are M different offset error signals and M different adaptive algorithms operating in conjunction with M different DACs providing offset control signals to M different ADCs.

The differences in the gain values of the ADCs produce an unwanted signal called the gain spurs. The frequencies of these spurs are ±F_(in)+kF₅/m where F_(in) represents a set of frequencies of the input signal. As can be seen from FIG. 1, the SNR variation is independent of the input signal frequencies F_(in). It, however, depends on the amplitude of the input signal. In order to obtain the gain error of each ADC, a certain window of length N_(g) samples is assumed. Each of the N_(g) samples from the output of each ADC are squared, and then a sum or average of these values is obtained from each ADC. Call each of the sums or averages as Y_(k) where k=1, 2, . . . , M. A single value representing the mean of these M values, say Y_(mean), is chosen as a reference gain value. The gain errors for the M different ADCs are obtained as Y_(k)−Y_(mean). The sign of each gain error, i.e., sign (Y_(k)−Y_(mean)), is used to drive an adaptive algorithm whose output represents a gain correction value for the corresponding ADC. The output from the adaptive algorithm is fed to a DAC whose output is a voltage or current that directly or indirectly controls the gain setting of each ADC. Thus, there are M different gain error signals and M different adaptive algorithms operating in conjunction with M different DACs providing gain control signals to M different ADCs.

The non-uniformity of the sampling instants of each ADC with respect to the system sampling instants of the MCTIADC gives rise to sampling or phase spurs. These spurs occur at the same frequencies as those due to the gain error. However, spurs due to the gain errors are orthogonal to those due to phase errors. In order to obtain the phase error, one first obtains the correlation between the samples of two adjacent ADCs. In other words, the samples of ADC₁ are correlated with the samples of ADC₂, samples of ADC₂ are correlated with the samples of ADC₃, and so on. The samples of ADC_(M) are correlated with samples of ADC₁ in the following cycle. These correlations are summed or averaged over a certain number of samples, N_(p). If Z_(k) denotes the sum or average of any correlation and Z_(mean) denotes the average of Z_(k), a phase error for any ADC can be formed as Z_(k), −Z_(mean). As in the case of offset and gain, Z_(mean) is assumed to be the reference phase value. Again, the sign of each phase error, i.e., sign (Z_(k)−Z_(mean)), is used to drive an adaptive algorithm whose output represents a phase correction value for the corresponding ADC. The output from the adaptive algorithm is fed to a DAC whose output is a voltage or current that directly or indirectly controls the phase setting of each ADC. Thus, there are M different phase error signals and M different adaptive algorithms operating in conjunction with M different DACs providing phase control signals to M different ADCs.

As mentioned above, the spur frequencies for gain and phase appear at ±F_(in)+kF₅/M. If one of the signal components of F_(in) is equal to

$\frac{{kF}_{s}}{2\; M},$

then it would be impossible to distinguish between that tone and the spur due to gain and phase mismatches. As a consequence of this, the algorithms for gain and phase correction tend to diverge. In order to circumvent this problem, a notch filter is introduced at the output of the ADC that will notch out these frequencies. The output from each notch filter is then used to evaluate the gain and phase errors, as mentioned above.

BRIEF DESCRIPTION OF THE DRAWINGS

The foregoing will be apparent from the following more particular description of example embodiments of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating embodiments of the present invention.

FIG. 1 illustrates SNR variation with input frequency of a typical Four-channel Time-Interleaved Analog-to-Digital Converter for various mismatch errors.

FIG. 2 is a schematic representing the M-channel Time-Interleaved Analog-to-Digital Converter (MCTIADC).

FIG. 3 is a spectrum of a single tone signal with offset mismatch error before correction in a four-channel time-interleaved ADC.

FIG. 4 is a schematic representing offset error calculation.

FIG. 5 is a schematic representing the recursive structure for effecting the offset adaptive algorithm.

FIG. 6 is a spectrum of a single tone signal with offset mismatch error after correction in a four-channel time-interleaved ADC.

FIG. 7 is a spectrum of a single tone signal with gain mismatch error before correction in a four-channel time-interleaved ADC.

FIG. 8 is a schematic representing gain error calculation.

FIG. 9 is a schematic representing the recursive structure for effecting the gain adaptive algorithm.

FIG. 10 is a spectrum of a single tone signal with gain mismatch error after correction in a four-channel time-interleaved ADC.

FIG. 11 is a spectrum of a single tone signal with phase mismatch error before correction in a four-channel time-interleaved ADC.

FIG. 12 is a schematic representing phase error calculation.

FIG. 13 is a schematic representing the recursive structure for effecting the phase adaptive algorithm.

FIG. 14 is a spectrum of a single tone signal with phase mismatch error after correction in a four-channel time-interleaved ADC.

FIG. 15 is a spectrum of a single tone signal with offset, gain and phase mismatch errors before correction in a four-channel time-interleaved ADC.

FIG. 16 is a spectrum of a single tone signal with offset, gain and phase mismatch errors after correction in a four-channel time-interleaved ADC.

FIG. 17 illustrates convergence of offset mismatch error in a four-channel time-interleaved ADC.

FIG. 18 illustrates convergence of gain mismatch error in a four-channel time-interleaved ADC.

FIG. 19 illustrates convergence of phase mismatch error in a four-channel is time-interleaved ADC.

FIG. 20 illustrates a spectrum of a multi-tone signal with offset, gain and phase mismatch errors before correction in a four-channel time-interleaved ADC.

FIG. 21 illustrates a spectrum of a multi-tone signal with offset, gain and phase mismatch errors after correction in a four-channel time-interleaved ADC.

FIG. 22 illustrates convergence of offset mismatch error in a four-channel time-interleaved ADC with a 100-tone signal.

FIG. 23 illustrates convergence of gain mismatch error in a four-channel time-interleaved ADC with a 100-tone signal.

FIG. 24 illustrates convergence of phase mismatch error in a four-channel time-interleaved ADC with a 100-tone signal.

FIG. 25 illustrates a magnitude response of the second-order notch filter.

FIGS. 26( a), 26(b) and 26(c) illustrate a multi-rate structure of a two-channel ADC with the notch filter.

FIG. 27 illustrates a magnitude response of the fourth-order notch filter.

FIGS. 28( a) and 28(b) illustrate a multi-rate structure of a four-channel ADC with the notch filter.

FIG. 29 illustrates a spectrum of a two-tone signal with a tone at Fs/8 before calibration.

FIG. 30 illustrates a spectrum of a two-tone signal with a tone at Fs/8 after calibration without the notch filter.

FIG. 31 illustrates a spectrum of a two-tone signal with a tone at Fs/8 after calibration with the notch filter.

FIG. 32 is an example communication device that may use the MCTIADC.

DETAILED DESCRIPTION OF THE INVENTION

A description of example embodiments of the invention follows.

A preferred embodiment deals with the estimation and correction of offset, gain and timing or phase mismatch errors in an M-channel Time-Interleaved Analog-to-Digital (MCTIADC) system. The estimation is done in the digital domain while the correction is performed in the analog domain. The various errors are estimated by performing signal processing operations on the output of all the ADCs while corresponding correction values are communicated to all the ADCs through Digital-to-Analog Converters (DACs). The output of each ADC is optionally passed through a notch filter to circumvent certain conditions in the input signal that will cause the algorithms to diverge. The details of the notch filter and its usefulness will be deferred to a later section below. The DACs provide appropriate voltages or currents and control either directly or indirectly the correction of each of the ADCs for the different mismatch errors.

FIG. 2 shows a schematic of an MCTIADC 200 wherein each of the M ADCs (210-1, 210-2, . . . , 210-M) is operating at a sampling rate of F₅/M and clocked at the appropriate respective phase φ_(k), for k=1 to M.

The commutator 230 operates at the sample rate F₅ and circles through the output of every ADC 210 to provide output y(n) at F₅. Outputs from each ADC 210 are input to a digital signal processor (DSP) 240 that performs the estimation of all the errors and provides analog outputs corresponding to offset, gain, and phase correction, represented by, O_(k), G_(k), and P_(k), respectively, to a corresponding one of the ADCs 210. Each of the ADCs may implement an offset, gain and phase correction as provided at the O_(k), G_(k), and P_(k), inputs, respectively, and may be implemented according to the charge domain pipeline ADCs described in the co-pending U.S. patent application Ser. No. 12/419,599, filed Apr. 7, 2009 already incorporated by reference above.

Also shown in FIG. 2 are the DACs 250 which convert the offset, gain and phase corrections O_(k), G_(k), and P_(k) to analog voltages before they are fed to the respective inputs of the ADCs 210. The DACs 250, of which there are 3×M in total (three DACs for each of the O, G, and P inputs to each of the M ADCs) can be implemented with any convenient circuit design to convert a digital input to an analog voltage, such as resistive ladder(s), look up table(s), amplifier(s), etc. that meets the desired clock rate and accuracy requirements. The DACs 250 and the optional notch filters 220 are described in more detail below.

Below is a detailed description of the estimation of offset, gain, and phase mismatch errors using the outputs of each ADC and their correction using adaptive algorithms that are performed within the DSP 240.

Offset Correction

Due to different offset values of the ADCs, offset spurs show up at kF₅/M frequencies. FIG. 3 shows the spectrum of a tone in a four-channel time-interleaved ADC sampling at 1 GHz where the offset spurs appear at fixed frequencies of 250 MHz and 500 MHz. In order to minimize the amplitude of these spurs, the offsets of each ADC must be determined. Towards this end, define

$\begin{matrix} {X_{k} = {\frac{1}{N_{o}}{\sum\limits_{n = 0}^{N_{o} - 1}\; {x_{k}(n)}}}} & (1) \end{matrix}$

where x_(k)(n) represents the samples from ADC_(k), N_(o) is the number of samples collected to obtain the average X_(k) and k=1, 2, . . . , M. Let

$\begin{matrix} {X_{mean} = {\frac{1}{M}{\sum\limits_{k = 1}^{M}\; X_{k}}}} & (2) \end{matrix}$

Now define an offset error for each ADC as

E _(k) ^(offset) =X _(k) −X _(mean)  (3)

for k=1, 2, . . . , M. It can be seen from the above equation that X_(mean) provides an overall reference value in each iteration so that an adaptive algorithm can be used to minimize E_(k) ^(offset).

Now it is possible to provide an adaptive algorithm to correct the offset error in each ADC based on E_(k) ^(offset), for k=1, 2, . . . , M.

Let ODAC_(k) be the one of the DACs 250 that provides the offset correction input O_(k) to ADC_(k). Let R_(o) be the size of the ODAC_(k). For example, for an 8-bit ODAC_(k), R_(o)=2⁸=256. A step size that controls the convergence of the adaptive algorithm is denoted by μ_(k) ^(i) for ADC_(k) at the ith iteration. The value of μ_(k) ^(i) is constrained to be in the range [μ_(k) ^(offsetmin),μ_(k) ^(offsetmax)]. Let O_(k) ^(i) be the jth value input to the ODAC_(k). For example, for an 8-bit ODAC_(k), the values of O_(k) ^(i) can vary between [−128,127] or between [0,255]. The constant O_(bias) is a value that allows the correction to be done with respect to a certain value. For instance, O_(bias)=R_(o)/2=128 when the input to the ODAC_(k) lies in the range [0,255]. On the other hand when the range of the ODAC_(k) input values is [−128,127], O_(bias) can assume a value of zero. Let α_(k) ^(i) denote a variable that provides correction to the ODAC_(k) input O_(k) ^(i) associated with ADC_(k) at the ith iteration. It is now possible to write the adaptive algorithm for offset correction as

$\begin{matrix} {O_{k}^{i} = {O_{bias} + {{round}\mspace{14mu} \left( \alpha_{k}^{i} \right)}}} & (4) \\ {\alpha_{k}^{i + 1} = {\alpha_{k}^{i} + {{sign}\mspace{14mu} \left( E_{k}^{offset} \right)\mu_{k}^{i}}}} & (5) \\ {\mu_{k}^{i + 1} = {{{\max\left( {\frac{\mu_{k}^{i}}{2},\mu_{k}^{{offset}\mspace{14mu} \min}} \right)}\mspace{14mu} {for}\mspace{14mu} i} = r_{k}}} & (6) \end{matrix}$

where α_(k) ⁰=0, μ_(k) ⁰=μ_(k) ^(offsetmax), and r_(k) is any arbitrary positive number. The convergence can be controlled by μ_(k) ^(i) by changing its value at every r_(k)th iteration.

A schematic that shows how the DSP 240 can perform the calculation of X_(k) (for k=1 to M) and X_(mean) is shown in FIG. 4. The output from each ADC_(k) is accumulated (with a corresponding summer 270-k and delay 272-k) for N_(o) samples to provide each X_(k). The port selector 260 then selects each accumulated ADC output X_(k) in turn, since the outputs from the accumulators are available in a time-interleaved fashion. The result is then further accumulated (by accumulator 280) and averaged (by multiplier 282), by dividing the accumulated sum by M to provide X_(mean). Next, X_(mean) is subtracted from each X_(k) to provide E_(k) ^(offset).

In FIG. 5, a schematic for an adaptive algorithm for how the DSP 240 can perform offset correction is depicted. The sign 310 of each E_(k) ^(offset) is first multiplied 315 by the adaptation step-size μ_(k) ^(i) and then accumulated 320. The accumulated value in each iteration is rounded 325 to the nearest integer value and added 330 to the offset bias, O_(bias), to provide the offset correction value O_(k) ^(i) to a corresponding one of the DACs, namely ODAC_(k). The output from ODAC_(k) directly or indirectly controls the offset setting on ADC_(k). Such an adaptive process converges to an optimal value that minimizes the offset error in each ADC.

FIG. 6 shows the spectrum of the tone mentioned in FIG. 3 after correction. As can be seen from the figure, the offset spurs at 250 MHz and 500 MHz are reduced.

Gain Correction

Gain differences in the ADCs produce gain spurs at ±F_(n)+kF_(s)/M frequencies, where F_(in) is the set of input frequencies and k=1, 2, . . . , M. FIG. 7 shows the spectrum of a 145 MHz tone in a four-channel time-interleaved ADC (without correction) sampling at 1 GHz where the gain spurs appear at 105 MHz, 355 MHz and 395 MHz. In order to reduce the amplitude of these spurs, the power of the signals from each ADC are determined. Towards this end, define

$\begin{matrix} {Y_{k} = {\frac{1}{N_{g}}{\sum\limits_{n = 0}^{N_{g} - 1}{x_{k}^{2}(n)}}}} & (7) \end{matrix}$

where X_(k) (n) represents the samples from ADC_(k), N_(g) is the number of samples collected to obtain Y_(k) and k=1, 2, . . . , M. Let

$\begin{matrix} {Y_{mean} = {\frac{1}{M}{\sum\limits_{k = 1}^{M}\; Y_{k}}}} & (8) \end{matrix}$

Now define a gain error for each ADC as

E _(k) ^(gain) =Y _(k) −Y _(mean)  (9)

for k=1, 2, . . . , M. It can be seen from the above equation that Y_(mean) provides a reference value for the power in each iteration so that an adaptive algorithm can be used to minimize E_(k) ^(gain). Below is outlined an adaptive algorithm to correct the gain error in each ADC based on E_(k) ^(gain), for k=1, 2, . . . , M.

Let GDAC_(k) be the one of the DACs 250 that provides the gain correction to ADC_(k). Let R_(G) be the size of the GDAC_(k). A step size that controls the convergence of the adaptive algorithm associated with gain correction is denoted by V_(k) ^(i) for ADC_(k) at the ith iteration. The value of v_(k) ^(i) lies in the range [v_(k) ^(offsetmin),v_(k) ^(offsetmax)]. Let G_(k) ^(i) be the value input to the GDAC_(k). Again, the values of G_(k) ^(i) can vary between [−128,127] or between [0,255] if R_(G)=256. The constant G_(bias) is a value that allows the correction to be done with respect to a certain value. For the case when G_(bias)=R_(G)/2=128, the input to the GDAC_(k) lies in the range [0,255]. On the other hand, when the range of the GDAC_(k) input values is in [−128,127], G_(bias)=0. Let β_(k) ^(i) denote a variable that provides correction to the GDAC_(k) input G_(k) ^(i) associated with ADC_(k) at the ith iteration. Now the adaptive algorithm for gain correction can be written as

$\begin{matrix} {G_{k}^{i} = {G_{bias} + {{round}\left( \beta_{k}^{i} \right)}}} & (10) \\ {\beta_{k}^{i + 1} = {\beta_{k}^{i} + {{{sign}\left( E_{k}^{gain} \right)}v_{k}^{i}}}} & (11) \\ {v_{k}^{i + 1} = {{{\max \left( {\frac{v_{k}^{i}}{2},v_{k}^{gainmin}} \right)}\mspace{14mu} {for}\mspace{14mu} i} = s_{k}}} & (12) \end{matrix}$

where β_(k) ⁰=0, v_(k) ⁰=v_(k) ^(gainmax), and s_(k) is an arbitrary positive number. The convergence can be controlled by V_(k) ^(i) by changing its value at every s_(k)th iteration.

A schematic that shows how the DSP can perform a calculation of Y_(k) and Y_(mean) is shown in FIG. 8. The output from each ADC_(k) is squared 810 and accumulated 820 for N_(g) samples to provide Y_(k). The port selector 830 then selects the squared and accumulated output of each ADC in turn and accumulates 840 the result. This is followed by an averaging operation 850 that is effected by dividing the accumulated sum by M to provide Y_(mean). Next, Y_(mean) is subtracted 860 from each Y_(k) to provide E_(k) ^(gain).

In FIG. 9, a schematic for an adaptive algorithm for the DSP 240 to perform gain correction is shown. The sign 910 of each E_(k) ^(gain) is multiplied 920 by the adaptation step-size and accumulated 930. The accumulated value in each iteration is rounded 940 to the nearest integer value and added 950 to the gain bias, G_(bias), to provide the gain correction value to GDAC_(k). The output from GDAC_(k) directly or indirectly controls the gain setting on GADC_(k).

The above adaptive process converges to an optimal value that minimizes the gain error in each ADC. FIG. 10 shows the spectrum of the tone mentioned in FIG. 7 after gain mismatch correction. As can be seen, the gain spurs at 105 MHz, 355 MHz and 395 MHz have been reduced.

Phase Correction

Since all the ADCs 210 do not have uniform sample instants in reference to the sampling frequency of the MCTIADC 200, timing or phase spurs show up at the same frequencies as those due to gain errors. One difference however is that gain spurs are orthogonal to the phase spurs. FIG. 11 shows the spectrum of a 145 MHz tone in an uncorrected four-channel time-interleaved ADC sampling at 1 GHz with phase spurs. As can be seen, the phase spurs occur at the same frequencies as those shown in FIG. 7. In order to minimize the amplitude of these spurs, a cross-correlation between any two adjacent ADCs must be determined. In view of this, define

$\begin{matrix} \begin{matrix} \left. {Z_{k} = {{\frac{1}{N_{P}}{\sum\limits_{n = 1}^{N_{p}}\; {x_{k}\left( {n - 1} \right)}}} - {x_{k + 1}\left( {n - 1} \right)}}} \right)^{2} & {{{{{for}\mspace{14mu} k} = 1},2,\ldots \mspace{14mu},{M - 1}}} \\ \left. {= {{\frac{1}{N_{P}}{\sum\limits_{n = 1}^{N_{p}}\; {x_{M}\left( {n - 1} \right)}}} - {x_{1}(n)}}} \right)^{2} & {{{{for}\mspace{14mu} k} = M}} \end{matrix} & (13) \end{matrix}$

where X_(k)(n) represents the samples from ADC_(k), N_(p) is the number of samples collected to obtain the average Z_(k) and k=1, 2, . . . , M. Let

$\begin{matrix} {Z_{mean} = {\frac{1}{M}{\sum\limits_{k = 1}^{M}\; Z_{k}}}} & (14) \end{matrix}$

Now define a phase error for ADC_(k) as

E _(k) ^(phase) =Z _(k) −Z _(mean)  (15)

for k=1, 2, . . . , M. In a manner similar to offset and gain error estimation, Z_(mean) provides a reference value for timing in each iteration so that an adaptive algorithm can be used to minimize E_(k) ^(phase). It is now possible to provide an adaptive algorithm to correct the phase error in each ADC based on E_(k) ^(phase), for k=1, 2, . . . , M.

Let PDAC_(k) be the DAC 250 that provides the timing or phase correction to ADC_(k). Let R_(p) be the size of the PDAC_(k). A step size that controls the convergence of the adaptive algorithm associated with phase correction is denoted by ξ_(k) ^(i) for ADC_(k) at the ith iteration. The value of ξ_(k) ^(i) is constrained to be in the range [ξ_(k) ^(phasemin),ξ_(k) ^(phasemax)]. Let P_(k) ^(i) be the value input to the PDAC_(k). Similar to the bias values in the offset and gain adaptive algorithms, the constant P_(bias) is a value that allows the correction to be done with respect to a certain value. Let γ_(k) ^(i) denote a variable that provides correction to the PDAC_(k) input P_(k) ^(i) associated with ADC_(k) at the ith iteration. It is possible to now write the adaptive algorithm for phase correction as

P_(k)^(i) = P_(bias) + round(γ_(k)^(i)) γ_(k)^(i + 1) = γ_(k)^(i) + sign(E_(k)^(gain))ξ_(k)^(i) $\xi_{k}^{i + 1} = {{{\max \left( {\frac{\xi_{k}^{i}}{2},\xi_{k}^{gainmin}} \right)}\mspace{14mu} {for}\mspace{14mu} i} = t_{k}}$

where γ_(k) ⁰=0, ξ_(k) ⁰=ξ_(k) ^(gainmax), and t_(k) is any arbitrary positive number. The convergence of the adaptive algorithm is controlled by ξ_(k) ^(i) by changing its value at every s_(k)th iteration.

A schematic that shows the calculation of Z_(k) and Z_(mean) by DSP 250 is shown in FIG. 12. For all ADC_(k), where k=1, 2, . . . , M−1, the squared 1210 difference 1220 of the delayed 1230 outputs from any two adjacent ADCs is accumulated for N_(p) samples to provide the corresponding X_(k). For ADC_(m), the present sample on ADC₁ is subtracted from the delayed input from ADC_(m) and squared. This value is then accumulated 1250 to obtain X_(M). The port selector 1260 then selects each ADC in turn and accumulates 1270 the result and performs an averaging operation by dividing 1280 the accumulated sum by M to provide Z_(mean). Next, Z_(mean) is subtracted 1290 from base each Z_(k) to provide E_(k) ^(phase).

In FIG. 13, a schematic for adaptive algorithm performing phase correction is shown. The sign 1310 of each E_(k) ^(phase) is multiplied 1320 by the adaptation step-size ξ_(i) ^(k) and accumulated 1330. The accumulated value in each iteration is rounded 1340 to the nearest integer value and added 1350 to the phase bias, P_(bias), to provide the phase correction value to PDAC_(k). The output from PDAC_(k) directly or indirectly controls the phase setting on ADC_(k).

The adaptive algorithm for phase correction converges in a manner that minimizes the magnitude of E_(k) ^(phase). FIG. 14 shows the spectrum of the tone mentioned in FIG. 3 after phase correction. As can be seen from the figure, the phase spurs at 105 MHz, 355 MHz and 395 MHz have been reduced.

So far what has been described are the adaptive algorithms pertaining to specific mismatch errors. In the presence of all the mismatches, viz., offset, gain and phase mismatches, the adaptive algorithms are run on a round-robin basis, starting with offset, gain and then phase. FIG. 15 shows the spectrum of a tone with all the mismatch errors while FIG. 16 shows the spectrum after mismatch errors have been reduced. As can be seen from the figure, the offset spurs at 250 MHz and 500 MHz, as well as gain and phase spurs at 105 MHz, 355 MHz and 395 MHz have been reduced.

FIGS. 17, 18, and 19 show the convergence of offset, gain and phase errors using the respective adaptive algorithms. It should be understood that the corrections can be made simultaneously, with one or more DSPs or hardware circuits, or by sharing a single DSP or a single hardware circuit in a round-robin fashion. What is important is that the expected rate of change in gain, offset and phase be slower than the rate at which the adaptive algorithms are performed.

The adaptive algorithms described thus far have shown to work for the case when the input is a single tone. It will be shown that the same set of algorithms will work for the case when the input signal is a wide-band signal. The difference between the two scenarios is the convergence time. The error functions tend to become non-linear and consequently the step size in the adaptive algorithms described would need to be much smaller than when the input signal is a single tone. FIG. 20 shows the spectrum of a wide-band signal comprised of many sinusoids in presence of offset, gain and phase mismatch errors while FIG. 21 shows the spectrum of the same signal after the spurs due to the mismatch errors have been reduced. FIGS. 22, 23, and 24 show the convergence of offset, gain and phase errors using the adaptive algorithms described above.

Notch Filter and its Usefulness

It was mentioned earlier that the spur frequencies for gain and phase appear at ±F_(in)+kF_(s)/M . If one of the signal components of F_(in) is equal to

$\frac{{kF}_{s}}{2M},$

then it is impossible to distinguish between this tone and the spurs due to gain and phase mismatches. Call such a frequency the gain-phase problem frequency (GPPF). As a consequence GPPF, the algorithms for phase and gain correction tend to diverge. In order to circumvent the divergence of these algorithms in the cases where the input spectrum has a GPPF, a notch filter is introduced at the output of the ADC that will notch out this frequency. The output from each notch filter can then used to evaluate offset, gain and phase errors in a way described above.

In order to develop the design of a notch filter to address this problem, first consider M=2. In this case, the GPPF is at Fs/4. Without any loss of generality, consider the design of a second-order notch filter that performs the notch filtering at these frequencies. A digital notch filter can be realized as

$\begin{matrix} {{G(z)} = {\frac{1}{2}\left( {1 + {A(z)}} \right)}} & (19) \end{matrix}$

where A(z) is an all-pass filter. The characteristics of G(z) are such that)

G(e ^(j0))=G(e ^(jπ))=1  (20))

G(e ^(jω) ⁰ )=0

where ω₀ is the angular notch frequency. A second-order transfer function to effect the all-pass filter is given by

$\begin{matrix} {{A(z)} = \frac{k_{2} + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + z^{- 2}}{1 + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + {k_{2}z^{- 2}}}} & (21) \end{matrix}$

where k and k₂ are multipliers defining the notch parameters. It can be shown that this choice of all-pass filter allows the independent tuning of ω₀ and the 3-dB bandwidth according to

$\begin{matrix} {k_{1} = {- {\cos \left( \omega_{0} \right)}}} & (22) \\ {k_{2} = \frac{1 - {\tan \left( {\Omega/2} \right)}}{1 + {\tan \left( {\Omega/2} \right)}}} & (23) \end{matrix}$

where Ω is the 3-dB bandwidth. Using Eqns. 19 and 21, one gets

$\begin{matrix} \begin{matrix} {{G(z)} = {\frac{1}{2}\left( {1 + \frac{k_{2} + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + z^{- 2}}{1 + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + {k_{2}z^{- 2}}}} \right)}} \\ {= {\frac{1 + k_{2}}{2}\frac{1 + {2k_{1}z^{- 1}} + z^{- 2}}{1 + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + {k_{2}z^{- 2}}}}} \\ {= {K\frac{1 + {2k_{1}z^{- 1}} + z^{- 2}}{1 + {{k_{1}\left( {1 + k_{2}} \right)}z^{- 1}} + {k_{2}z^{- 2}}}}} \end{matrix} & (24) \end{matrix}$

where K=(1+k₂)/2 is a scaling factor based on the value of k₂. First consider the case of a notch frequency at Fs/4, i.e. ω₀=762. It can be seen from Eqn. 22, that multiplier k₁=0. Having eliminated the need for a multiplier, Eqn. 24 can be written as

$\begin{matrix} \begin{matrix} {{G(z)} = {K\frac{1 + z^{- 2}}{1 + {k_{2}z^{- 2}}}}} \\ {= {H\left( z^{2} \right)}} \end{matrix} & (25) \end{matrix}$

where

$\begin{matrix} {{H(z)} = {K\frac{1 + z^{- 1}}{1 + {k_{2}z^{- 1}}}}} & (26) \end{matrix}$

From Eqn. 26, the impulse response can be written as

$\begin{matrix} {{h(n)} = {{\frac{K}{k_{2}}{\delta (n)}} + {\frac{K\left( {k_{2} - 1} \right)}{k_{2}}\left( {- k_{2}} \right)^{n}{u(n)}}}} & (27) \end{matrix}$

Now write Eqn. 23 in terms of Ω. After some manipulation, one gets

$\begin{matrix} {\Omega = {2{\tan^{- 1}\left( \frac{1 - k_{2}}{1 + k_{2}} \right)}}} & (28) \end{matrix}$

As can be seen from the above equation, the bandwidth depends on the value of k₂. Thus by appropriately choosing the value of k₂ as a Canonic Signed Digit (CSD) number, the need for a multiplication can be circumvented. For example, by choosing k₂=1−2⁻³=0.875, one can obtain a bandwidth of Ω=0.0424π. It must be noted that as k₂ approaches unity, K approaches unity. Hence, in many applications, the scale factor K can also be eliminated. Using K=1, the frequency response of a notch filter G(z) with the above value of k₂ is shown in FIG. 25 for ω₀=π/2.

For M=2, consider the two-channel interleaved ADC in conjunction with the notch filter characterized by H(z²) from a multi-rate signal processing point of view. FIG. 26( a) shows such a structure where equivalent output from each ADC is the output from a 2× decimator 2610. The commutator action is represented by the 2× interpolators 2620, the delay element 2630 and the adder 2632. The signals at the output of the adder 2632 and the notch filter 2635 operate at F_(s). By shifting the filter 2635 before the adder 2632, one gets the structure shown in FIG. 26( b). Finally, by using a property called Nobel Identity in multi-rate signal processing, one can move the filter 2635, given by H(z), before the interpolators 2620. Hence, starting with a second-order notch filter with a notch at ω₀=π/2, in FIG. 26( c) it has now been converted to a single pole filter running at Fs/2 rate. It is interesting to note that H(z) now has a notch at the Nyquist frequency of each ADC.

Now consider the case of M=4. The GPPF are at Fs/8 and 3F_(s)/8. In terms of normalized frequency these frequencies are π/4 and 3π/4. Consider an 2× upsampled version of G(z) given by Eqn. 25. One gets

$\begin{matrix} \begin{matrix} {{G\left( z^{2} \right)} = {K\frac{1 + z^{- 4}}{1 + {k_{2}z^{- 4}}}}} \\ {= {H\left( z^{4} \right)}} \end{matrix} & (29) \end{matrix}$

The magnitude response of this filter, viz., G(z²) is shown in FIG. 27. It is evident from the magnitude response that the notch frequencies are π/4 and 3π/4. Again, let us look at a four-channel ADC from a multi-rate perspective with the notch filter given by H(z⁴). FIG. 28( a) shows the four-channel structure while FIG. 28( b) shows the equivalent structure modified by the property of the Noble Identity. Again, it can be seen the H(z) is the basic filter that provides the notch filtering of all the GPPF.

Extending this theory to the case of M-channel interleaved ADC, one can arrive at the structure shown in FIG. 2 wherein each notch filter characterized by Eqn 26 is used.

A simulation demonstrated the usefulness of the notch filter in a four-channel time-interleaved ADC. In this simulation, a two-tone signal was considered, with one tone at 125 MHz, i.e. at Fs/8 where F_(s)=1 GHz, and the other at any arbitrary location. Here the other tone was chosen at 45.123 MHz. FIG. 29 shows the spectrum of the signal with offset, gain and phase spurs. As can be seen from that figure, for the input tone at F_(s)/8, the gain and phase spur appear at F₅/8, F_(s)/4 and 3F_(s)/8. The input tone at 45.123 MHz produces gain and phase spurs at approximately 205 MHz, 295 MHz and 455 MHz. There is also a tone at 250 MHz due to offset spur. FIG. 30 shows the spectrum after 15000 iterations of the algorithm. Here the notch filter was not used. As can be seen from the figure, the offset algorithm performs fairly well and hence the tone at 250 MHz has been reduced. However, all the other spurs still exist.

FIG. 31 shows the spectrum after 15000 iterations where the notch filter has been used. As can be seen from the figure, all the spurs have been reduced significantly.

High sample rate, time interleaved ADCs such as that described above can find application in many different types of systems. One such application is the receivers used in communication systems. Such receivers have historically used analog tuner devices to demodulate a small portion of the input signal spectrum down to a low frequency. Relatively speaking, the tuner output has a low center frequency and low total bandwidth, thus allowing a low speed analog-to-digital converter to be used to digitize the data.

Certain popular communication system receivers such as that used in cable modem and set-top-box systems are trending toward processing more channels to provide faster broadband access and more video services to the home. One alternative to having a tuner for each channel is to digitize the entire bandwidth of the cable system. Once this is achieved, the number of channels decoded from the spectrum is completely defined in the digital domain. Thus the incremental cost for each additional channel is relatively low and should decrease rapidly over time as digital process technology advances. Digitizing the entire bandwidth of the cable system requires a very high sampling rate; therefore, an interleaved system may provide advantages over other conversion techniques. The increasing need for a wideband spectrum incorporating frequency multiplexed signals makes cable and other communication systems an excellent application of this invention.

FIG. 32 shows an example communication device, such as cable gateway 3100 connected to a cable network 3108, which may be a coaxial, optical fiber, or hybrid fiber/coaxial cable television (CATV) network. The cable gateway 3100 transmits data to and receives data from customer premises equipment 3112. Typically, customer premises equipment 112 includes computers, televisions, and telephones. The cable gateway 100 disclosed herein can be configured to operate according to any suitable specification for transmitting and receiving data, including but not limited to DOCSIS 3.0, Comcast RNG, SCTE 40, T3/S10 ATSC, or OpenCable specifications. Certain specifications require cable modems and cable gateways to tune multiple channels at the same time for receiving television, voice, and data signals. (For example, DOCSIS 3.0 specifies the ability to independently tune at least four channels.) The ability to tune multiple channels is also necessary to watch different television channels on different televisions.

Unlike conventional cable gateways, the cable gateway 3100 shown in FIG. 32 uses a wideband, multi-channel, time-interleaved, analog-to-digital converter (MCTIADC) 3206 to digitize signals received from the cable network 3108. The output from the wideband ADC 3206 can be tuned digitally, rather than with analog tuners, resulting in lower power consumption compared to alternative methods. More particularly, in the example cable gateway 3100, signals transmitted to and from the cable network 3108 are coupled via a diplexer 3202, which separates downstream signals 3220 from upstream signals 3222. In general, CATV networks are asymmetric networks: the bandwidth dedicated to the downstream signals 3220 is greater than the bandwidth dedicated to the upstream signals 3222.

The diplexer 3202 directs downstream traffic to a variable-gain amplifier (VGA) 3204, which amplifies the received signal before transmitting it through a filter 3205 to a wideband time-interleaved ADC 3206. The time-interleaved ADC 3206 digitizes the received signal, then passes the digitized downstream signals 3240 to a digital tuner and quadrature-amplitude-modulation (QAM) or other type demodulator 3208. (Alternative embodiments may use other suitable modulation schemes.) In some embodiments, the digital tuner and QAM demodulator 3208 tunes and demodulates the amplified, filtered, and digitized downstream signals 3240 in accordance with either 64-QAM or 256-QAM techniques to recover the underlying information.

While this invention has been particularly shown and described with references to example embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the scope of the invention encompassed by the appended claims. 

1. An apparatus comprising: a clock signal generator, for generating a plurality, M, of clock signals at a frequency f and a period T; each of the M clock signals having a different one of a selected plurality, M, of clock phases a plurality, M, of Analog to Digital Converters (ADCs) coupled to the clock signal generator, the ADCs for converting an input signal to a set of ADC outputs as M digital values in response to a respective one of the M clock signals, each of the ADCs having an offset correction input, a gain correction input, and a phase correction input; an adaptive processor, coupled to receive the M digital values, the adaptive processor estimating one or more correction signals for at least one of offset, gain, and phase error in at least one of the MADCs, the adaptive processor determining the correction signals by: determining M accumulated values, X_(k), for k=1 to M, by individual accumulation of the M digital values over a predetermined number of ADC output samples; determining a reference value, X_(mean), from a combination of the M accumulated values; providing an adjusted set of set of digital values, E^(offset)k, for k=1 to M, from the M accumulated values, X_(k), and the reference value X_(mean), from the adjusted set of digital values, determining at least one of an offset, gain and phase correction value corresponding to one or more estimated correction signals to be applied to correct at least one of offset, gain, and phase error of at least one of the ADCs; the estimated correction signals connected to at least one of the offset, gain, and/or correction inputs of the ADCs; and a multiplexer, for interleaving the M digital values output by the ADCs to form a digital representation of the input signal.
 2. The apparatus of claim 1 additionally comprising: one or more notch filters, coupled to the ADCs, to filter spur frequency content in the M digital values resulting from at least one of the gain or phase error.
 3. The apparatus of claim 1 wherein the adaptive processor further determines the offset correction values by determining an accumulated offset for each of the plurality, M, of ADCs.
 4. The apparatus of claim 1 wherein the adaptive processor further determines the gain correction value by determining an accumulated gain error for each of the plurality, M., of ADCs.
 5. The apparatus of claim 1 wherein the adaptive processor further determines the phase correction value by determining an accumulated phase error for each of the plurality of ADCs.
 6. The apparatus of claim 1 wherein the adaptive processor sequentially determines the offset, gain, and phase correction values.
 7. The apparatus of claim 1 further including a plurality of digital-to-analog converters (DACs), where each DAC converts at least one of an offset, gain, or phase correction value to an analog signal input to at least one of the offset, gain, or phase inputs to the plurality, M, of ADCs.
 8. The apparatus of claim 1 further including additional adaptive processors, the adaptive processors configured to process offset, gain, and phase correction values in parallel.
 9. The apparatus of claim 1 wherein the adaptive processor further determines the offset correction value and further comprising: a signum block for determining a sign of the offset error; a multiplier for multiplying the output of the signum block by an step size; a feedback loop for summing and delaying an output of the multiplier; a rounding block for rounding an output of the feedback loop; and an offset biasing block for biasing the output of the rounding block by an offset bias amount.
 10. The apparatus of claim 1 wherein the adaptive processor further determines the gain correction value and further comprises: a signum block for determining a sign of the gain error; a multiplier for multiplying the output of the signum block by an gain step size; a feedback loop for summing and delaying an output of the multiplier; a rounding block for rounding an output of the feedback loop; and a gain biasing block for biasing the output of the rounding block by a gain bias amount.
 11. The apparatus of claim 1 wherein the adaptive processor further determines the phase correction value and further comprises: a signum block for determining a sign of the phase error; a multiplier for multiplying the output of the signum block by a phase step size; a feedback loop for summing and delaying an output of the multiplier; a rounding block for rounding an output of the feedback loop; and a phase biasing block for biasing the output of the phase rounding block by a phase bias amount.
 12. The apparatus of claim 1 additionally comprising: a receiver, for receiving the input signal from a communication system.
 13. The apparatus of claim 1 additionally comprising: a diplexer, coupled to receive a cable signal from a cable network, and to provide a received cable signal; a variable gain amplifier, for amplifying the received cable signal and to provide an amplified cable signal; and a filter, connected to filter the amplified cable signal to provide the input signal to the plurality, M, of ADCs.
 14. A method comprising the steps of: generating a plurality, M, of clock signals at a frequency f and a period T; each of the M clock signals having a different one of a selected plurality, M, of clock phases; converting an input analog signal to M digital signals via M individual ADC operations, to generate a set of ADC outputs as M digital values each in response to a respective one of the M clock signals, each of the M ADC operations having an offset correction input, a gain correction input, and a phase correction input; estimating one or more correction signals for at least one of offset, gain, and phase error in at least one of the M ADC operations, via the further steps of: determining M accumulated values, Xk, for k=1 to M, by individual accumulation of the M digital values over a predetermined number of ADC operation output samples; determining a reference value, Xmean, from a combination of the M accumulated values; providing an adjusted set of set of digital values, Eoffsetk, for k=1 to M, from the M accumulated values, Xk, and the reference value Xmean; from the adjusted set of digital values, determining at least one of an offset, gain and phase correction value corresponding to one or more estimated correction signals to be applied to correct at least one of offset, gain, and phase error of at least one of the ADC operations; coupling the estimated correction signals to at least one of the offset, gain, and/or correction inputs of the individual ADC operations; and interleaving the M digital values output by the ADCs to form a digital representation of the input signal.
 15. The method of claim 14 additionally comprising: applying one or more notch filter operations to the outputs of the ADCs, to filter spur frequency content in the M digital values resulting from at least one of the gain or phase error.
 16. The method of claim 14 further comprising: determining the offset correction values by determining an accumulated offset for each of the plurality, M, of ADCs.
 17. The method of claim 14 further comprising: determining the gain correction value by determining an accumulated gain error for each of the plurality of ADCs.
 18. The method of claim 14 further comprising: determining the phase correction value by determining an accumulated phase error for each of the plurality of ADCs.
 19. The method of claim 14 further comprising; sequentially determining the offset, gain, and phase correction values.
 20. The method of claim 14 further comprising: Digital to Analog converting at least one of an offset, gain, or phase correction value to one or more analog input signals; and coupling the analog input signals to at least one of the offset, gain, or phase inputs of the plurality, M, of ADCs.
 21. The method of claim 14 further comprising: processing the offset, gain, and phase correction values in parallel.
 22. The method of claim 14 further comprising: determining a sign of the offset error; multiplying the output of the signum block by an step size; summing and delaying an output of the multiplier; rounding an output of the feedback loop; and biasing the output of the rounding block by an offset bias amount.
 23. The method of claim 14 further comprising: determining a sign of the gain error; multiplying the output of the signum block by an gain step size; summing and delaying an output of the multiplier; rounding an output of the feedback loop; and biasing the output of the rounding block by a gain bias amount.
 24. The method of claim 14 further comprising: determining a sign of the phase error; multiplying the output of the signum block by a phase step size; summing and delaying an output of the multiplier; rounding an output of the feedback loop; and biasing the output of the phase rounding block by a phase bias amount.
 25. The method of claim 14 additionally comprising: receiving the input signal from a communication system.
 26. The method of claim 14 additionally comprising: receiving a cable signal from a cable network, to provide a received cable signal; amplifying the received cable signal to provide an amplified cable signal; and filtering the amplified cable signal to provide the input signal to the plurality, M, of ADCs. 